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Abstract 

Causality alone suffices to set a lower bound on the period of ro- 
tation of relativistic stars as a function of their maximum observed 
mass. That is, by assuming a one-parameter equation of state (EOS) 
that satisfies v SO und < c an d that allows stars with masses as large 
as the largest observed neutron-star mass, M™% x , we find P[ms] > 
0.282 + 0.196 (M™l x /M & - 1.442V The limit does not assume that 



the EOS agrees with a known low-density form for ordinary matter, 
but if one adds that assumption, the minimum period is raised by a 
few percent. Thus the current minimum period of uniformly rotat- 
ing stars, set by causality, is 0.28ms (0.29ms for stars with normal 
crust). The minimizing EOS yields models with a maximally soft 
exterior supported by a maximally stiff core. An analogous upper 
limit set by causality on the maximum mass of rotating neutron stars 
requires a low-density match and the limit depends on the match- 
ing density, e m . We recompute it, obtaining a slightly revised value, 
M™? x ~ 6.1 (2 x 10 14 g/cm 3 /e m ) 1/2 M . 
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1 INTRODUCTION 



The upper limit set by gravity on the rotation of neutron stars is sensitive 
to the equation of state (EOS) above nuclear density, and the current large 
uncertainty in the EOS implies a large uncertainty in the maximum rate 
of rotation. One can, however, find an upper limit on rotation which is 
independent of the EOS. This limit is set by causality together with the 
requirement that the EOS allows stars with masses as large as the largest 
observed neutron-star mass. 

Glendenning (1992) first estimated this causally limited period for gravi- 
tationally bound stars. He considered a flexible two-parameter, gamma-law 
ansatz for the EOS at high energy density; by varying the two parameters, 
he found equations of state (EOSs) which appeared to minimize the period 
as a function of the maximum nonrotating (or spherical) mass M™% x allowed 
by the EOS. Our work improves upon his estimate in several ways. First, 
where Glendenning used an empirical formula to estimate the minimum pe- 
riod from a set of nonrotating models, we construct models using a code 
that constructs rapidly rotating relativistic stars. Second, our numerical in- 
vestigation singles out a simple analytic form for the EOS which minimizes 
the period for a given M™ff. This minimum-period EOS has a simple phys- 
ical interpretation; it allows the most centrally condensed star while still 
supporting a mass M™% x . Specifically, we find that if M™% x = 1.442 M 
(currently the largest accurately measured mass of neutron stars (Taylor & 



3 



Weisberg 1989), the minimum period allowed is 0.28ms. This minimum pe- 
riod is about 13 % less than that estimated by Glendenning (1992). Third, 
we emphasize that there exists an upper limit on rotation set only by causal- 
ity, which requires no matching to a low density EOS. Finally, we find an 
exact scaling between the minimum period P min and M™™. The scaling is 
due to the specific form of the minimum-period EOS, and implies that the 
central energy densities, masses, and radii of the maximum-mass spherical 
(and maximum-mass rotating) stars scale between different minimum-period 
EOSs. 

The maximum mass of gravitationally bound stars set by causality is, in 
contrast, quite sensitive to the matching density (Hartle & Sabbadini 1977) 
(the initial Rhoades and Ruflini (1974) limit was computed for a particu- 
lar choice of matching density). For rotating neutron stars, the maximum 
mass set by causality was first obtained by Friedman and Ipser (1987). This 
maximum mass has not, however, however been recomputed since recent im- 
provements in relativistic, rotating star codes, and we do so here. Using 
the FPS EOS (Pandharipande & Ravenhall 1989) at low density, we find 
M max = 6.1M (e m /2 x 10 14 g /cm 3 ) -1 / 2 , where e m is the matching energy 
density. 
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2 ASSUPMTIONS FOR THE EQUATION 
OF STATE 



For number densities n < 0.1 fm 



-3 



(slightly less than nuclear saturation 



density) the FPS EOS (Pandharipande & Ravenhall 1989) is expected to 
describe normal matter with reasonable accuracy (Ravenhall 1995). Above 
this density, however, current observations are compatible with a variety of 
proposed equations of state that span a large range of compressibilities and 
predict correspondingly different minimum periods. One can, however, use 
causality to set a lower limit on the rotation period of stable relativistic stars 
that is independent of the high density EOS. More precisely, we assume: 

(0) A relativistic star is described as a self- gravitating, uniformly rotating 
perfect fluid with a one-parameter EOS. 

Cold neutron stars satisfy this assumption to high accuracy (Friedman 



Here v soun d is the phase velocity of sound waves, p the pressure, and 
e the energy density of the stellar perfect fluid matter. Relativistic 
fluids are governed by hyperbolic equations whose characteristics lie 
outside the light cone (and thus their initial value formulation violates 
causality) if v soun d > 1 (here and throughout the paper we have set 
c = 1) (Geroch & Lindblom 1991). That is, if both the equilibrium 



& Ipser 1992). 
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star and its small oscillations are described by a one-parameter EOS, 
then assumption (1) is implied by causality. 

These two assumptions alone yield an upper limit on rotation as a function 
of the maximum observed mass of a relativistic star. Only a slightly more 
stringent limit is obtained if one makes the further asumption: 

(2) The EOS is known (matches FPS) at low density. 

Relativistic stars with normal crusts are thought to be accurately de- 
scribed by the FPS EOS up to number densities n « 0.1 fm~ 3 . We 
adopt n m = 0.1 fm~ 3 as a conservative value for the number density 
n m at which the EOS at high density must match the FPS EOS for 
low densities. 

Below we examine limits on the rotation of stable equilibrium models cor- 
responding first to assumptions (0) and (1) only, and then to assumptions 
(0)-(2); as noted, we find that the additional assumption of a match to the 
known low-density equation of state has a small effect on the minimum period 
computed using (0) and (1) only. 

3 EQUATION OF STATE YIELDING THE 
UPPER LIMIT ON ROTATION 

A uniformly rotating, gravitationally bound star rotates with a period greater 
than or equal to its Kepler period. The Kepler period is the period of a free 
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particle in circular orbit at the star's equator. A soft EOS yields stellar 
models with dense central cores, large binding energies, and thus smaller 
Kepler periods than models built from stiff EOSs. Soft EOSs, however, 
cannot support massive stars. This suggests that models with minimum 
periods arise from EOSs which are stiff at high density, allowing stiff cores to 
support against collapse, but soft at low density, allowing small radii and thus 
fast rotation. Our numerical results are consistent with this expectation, and 
strongly suggest a simple form for the EOS yielding minimum period stars 
and satisfying the assumptions listed in section |2|. We distinguish the two 
instances where the EOS does or does not match smoothly to the FPS EOS 
at low density and satisfy assumption (2) above. 

3.1 Minimum-Period EOS 

Considering only assumptions (0) and (1), so that the EOS is constrained 
only by causality, the EOS yielding the minimum period stars for a given 
maximum spherical star mass M™^ x is divided into two regions. Above a 
density ec the EOS is maximally stiff with dp/de = 1, while for e < ec the 
EOS is maximally soft with vanishing pressure p = 0: 



This EOS is depicted in Figure |TJ. It gives a stringent upper limit set by 
causality on the rotation of uniformly rotating, gravitationally bound stars, 
independent of any specific knowledge about the EOS for the matter com- 




(1) 
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posing the star. The energy density ec is the single free parameter of this 
minimum-period EOS. By choosing ec one selects the maximum spherical 
star mass M™fi^ that the minimum-period EOS allows. Further, a one to one 
correspondence exists between ec and other characteristics of the minimum- 
period EOS, such as the maximum mass of rotating stars. 

3.2 Minimum-Period EOS with FPS Low Density EOS 

Considering assumptions (0), (1) and (2) above, so that the EOS matches 
the FPS EOS at low density, the minimum-period EOS is divided into three 
regions. Above a density ec the minimum-period EOS is maximally stiff with 
dp/de = 1. Between the matching density to the FPS EOS e m and ec the 
EOS is maximally soft with dp/de = 0. Below e m the EOS is given by the 
FPS EOS. Explicitly, 

' PFPs(e) e < e m , 

p(e) = < p m e m <e< e c , (2) 

p m + e-ec e> e C - 
This EOS is also depicted in Figure [l]. The maximally soft, constant pressure 
region corresponds to a first-order phase transition for a single component 
system. For a fixed low density EOS and matching energy density e m , the 
causal limit energy density ec parameterizes the class of minimum-period 
EOSs. Again, by choosing ec one selects the maximum spherical star mass 
M^f that the minimum-period EOS allows, and a one to one correspondence 
exists between ec and other characteristics of the minimum-period EOS. 
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Figure Q shows the minimum-period EOS (^) yielding a maximum mass 
spherical star with M™^ = 1.442 M . The energy density at which the EOS 
changes from zero pressure to dp/de = 1 is = 2.156 x 10 15 g/cm 3 , and 
the central energy density of the minimum-period star is 4.778 x 10 15 g/cm 3 . 
Also shown in Figure |] is the minimum-period EOS @ which matches to 
the FPS EOS at n m = 0.1 fm~ 3 and yields a maximum mass spherical star 
with M™% x = 1.442 M . For this EOS e c = 2.157 x 10 15 g/cm 3 and the 
central energy density of the minimum-period star is 5.274 x 10 15 g/cm 3 . The 
numerical work leading to the EOSs (|l|) and (|2]) is outlined in section |j. 

Figure Q shows the minimum period P m in of gravitationally bound stars 
as a function of the maximum allowed mass of nonrotating stars M™% x . The 
lower curve (which is a straight line) corresponds to the minimum-period EOS 
([!]) satisfying assumptions (0) and (1), while the upper curve corresponds to 
the minimum-period EOS (0) matched to the FPS EOS at n m = 0.1 fm~ 3 . 

The minimum-period EOS ([!]) yields stars entirely at the causal limit 
with a nonzero surface energy density ec- This surface energy density ec is 
the only free parameter and the only dimensionful parameter. The family of 
maximally rotating equilibrium stars yielded by the minimum-period EOS 
([]]) with different ec are characterized by this one dimensionful parameter; it 
follows that all properties of the maximally rotating stars scale according to 
their dimensions in gravitational units (with c = G = 1), [e] = [M~ 2 ], [P] = 
[M] = [R]. Thus, the following relations hold between different maximally 
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rotating stars computed from minimum-period EOSs with different ec- 
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where R is the equatorial radius. Numerically, we find the linear relation 
PminM = 0.200 M™r/M , or 

P ( M max 

= 0.288 + 0.200 - 1.442 

ms \ M Q 

(Note, however, that we have only considered the period of the maximum- 
mass rotating model. In general the maximum-mass model is distinct from 
the minimum-period model, but the difference in the period between the two 
models is small. See section ||.) 

The upper curve in Figure ^| corresponds to EOSs (H) which match the 
FPS EOS at low density. This curve is almost linear, since the matching 
number density n m = 0.1 fm~ 3 is low enough that the causal limit region 
of the EOS dominates the bulk properties of the star. That is, a minimum- 
period star with a low matching density mimics a star entirely at the causal 
limit and having a nonzero surface energy equal to ec- Hence the scaling 
relations above are nearly exact for minimum-period stars with low matching 
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density. Numerically we find for the EOSs (0) 

P ( M max 

= 0.295 + 0.203 - 1.442 

ms \ Mq 

which is linear to an accuracy better than 0.5 %. For M^ x = 1.442 M 
the minimum period is 0.29 ms. The above formulae allow one to calculate 
without intensive numerical computations the absolute minimum period if 
and when new observations revise the current M™% x = 1.442M Q limit. 

4 NUMERICAL EVIDENCE IN SUPPORT 
OF THE MINIMUM-PERIOD EOS 

As noted above, our numerical results are consistent with the expectation 
that the EOSs yielding the minimum-period stars are maximally stiff at high 
density and maximally soft at low density. Below we present the numerical 
evidence supporting our claim that the EOSs ([I]) and (Q) yield the maximum 
rotating stars, and detail how we searched the space of EOSs to find the 
minimum-period EOS 

4.1 Perturbations to the Minimum-Period EOS 

Evidence that the minimum-period EOS (0) yields the fastest rotating stars 
is shown in Figure |4]. Each point in Figure ^ represents an EOS obtained 
by making a "small" perturbation to the minimum-period EOS (||), and 
lies above the solid curve obtained from the minimum-period EOS. For any 
particular M™% x the perturbed EOSs yield maximally rotating stars which 
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rotate slower than the maximally rotating stars yielded by the minimum- 
period EOSs. 

We considered three different types of perturbations to the minimum 
period EOS. The first type of perturbed EOS considered is obtained from the 
minimum-period EOS by adding a "step function" to the minimum-period 
EOS : 

' PFPs(e) e < e m> 



P(e) 



Pr. 



p m + e-e 1 



p m + e 2 - ei 



e m < e < ei, 
ei < e < £2, 
£2 < e < ecr, 



(10) 



Pm + e 2 - ei + e - e c e > e c . 

By varying ei and 62, as well as ec, we sampled a large number of perturbed 
minimum-period EOSs. For e 2 — ei small, the minimum period for a given 
M™% x allowed by the perturbed EOS is slightly larger than that yielded by 
the minimum-period EOS (0). The minimum periods obtained using the 
"step function" perturbations are indicated by the open squares in Figure £|. 
For €2 — €1 large, the minimum periods yielded by the perturbed EOS are 
much larger and are not shown in Figure ^. Thus a small perturbation away 
from the minimum-period EOS @ increases the minimum period for a given 

]\/fmax 
IV1 sph ■ 

The second type of perturbation to the minimum-period EOS we consid- 
ered is obtained by decreasing the slope dp/ale away from unity for the causal 
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( e > £c) P ar t °f the minimum-period EOS. The filled triangles in Figure [| 
indicate the minimum periods obtained by perturbing the minimum-period 
EOS in this way. Again, the minimum period increases when the minimum- 
period EOS (^) is perturbed. The third type of perturbation we considered 
is obtained by increasing the slope dp/de away from zero for the constant 
pressure or e m < e < tc part of the minimum-period EOS. For e > ec the 
perturbed EOS again is at the causal limit with dp/de = 1. The crosses in 
Figure [| indicate the minimum periods obtained by perturbing the minimum- 
period EOS in this way, and again the minimum period increases when the 
minimum period EOS (Q) is perturbed. 

We have sampled a large number of perturbed minimum-period EOSs 
which satisfy assumptions (0),(1), and (2) of section |2|, and have found that 
the minimum period allowed by a perturbed EOS for a given M™% x is always 
greater than the minimum period yielded by the minimum-period EOS (0). 
These results, taken together with results obtained from a two parameter 
ansatz EOS described in the next section, strongly suggest that the minimum- 
period EOS (H) is the EOS which provides the upper limit set by causality 
for the rotation of gravitationally bound stars. 

4.2 Ansatz for Searching the Space of Equations of 
State 

To search the restricted space of EOSs which satisfy the assumptions in 
section ^ and to determine the class of EOSs which yield the minimum 
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period stars, we follow Glendenning (1992) and adopt a simple ansatz for 
the EOS. For low number density, the energy density and pressure are 

e(n) = e F p S , (11) 
p(n) = p F ps, (12) 

where €fps and pfps correspond to the FPS EOS. Above some number density 
n m , where the FPS EOS is suspect, the energy density and pressure are 

e(n) = p 7 
p(n) = Pr 

where e m and p 

density n m , which is typically 0.1 fm~ 3 . The two dimensionless parameters 
K and 7 parameterize the restricted EOS space. Assumptions (0) and (1) of 
section ^| require 

7 > 1, (15) 
k> 0, (16) 

k 7 < £ + 1. (17) 

Given k and 7, we use the parameterized form ( [TBI ) and ( |HD for densities 
greater than the matching density n m , but only if the EOS remains causal. 
If above some density nc the equation of state reaches the causal limit, so 
that $ > 1, then for n > nc the parameterized EOS is matched to an EOS 
which is always at the causal limit. In this instance the energy density and 
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Pin) = i |p c - e c + (Po + eo) (^) | , (19) 

where ec and are the energy density (|IciD and pressure ( fH] ) evaluated 
at n^. Note that the parameterization of pressure as a function of number 
density is arbitrary, since assures that $ = 1. If k and 7 are chosen 
so that the EOS remains causal then the parameterized form is used for 
all densities n > n m . Combining the different expressions for the different 

ranges of number density the EOS is 
cfpsH 
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p(n) 



n < n r . 



Pm { K [{^y] + 1 } n m <n<n c , (2 1) 

||pc-e c + (pc + ec) (^) 2 | n > n c , 
where the n > nc case applies only if nc exists and is finite. Note that the 
EOS is continuous, although its first derivative may be discontinuous at n m . 
We do not include a separate expression for the possibility of a p = constant, 
one-component first-order phase transition region, since this may be attained 



in the limit as k — > 0. One may search the restricted EOS space using (20) 
and (|2lD by varying the parameters k and 7. 
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The two-parameter EOS ( PUD and ( |2~T]) does not span the entire restricted 
EOS space. Still, as Figure [| shows, it is quite flexible, and spans a represen- 
tative subspace which includes both soft and stiff EOSs. This two-parameter 
family of EOSs proved sufficient to identify the class of EOSs yielding the 
minimum-period stars. 

4.3 Searching for the Minimum- Period Star 

We used a fully relativistic, numerical code to compute rapidly rotating mod- 
els of compact stars and find the maximum-mass model for a given EOS. Ster- 
gioulas (Stergioulas & Friedman 1995) wrote the code following the Cook, 
Shapiro and Teukolsky (1994) implementation of the Komatsu, Eriguchi and 
Hachisu (1989) method. Koranda made changes which substantially im- 
proved the speed of the code. 

Using the flexible ansatz (|20|) and (f21~f), we have searched for the class of 
EOSs which yield the minimum-period stars. As mentioned in the introduc- 
tion, we are interested in the minimum period as a function of M™^ x > the 
mass of the maximum mass nonrotating (spherical) star. We first used (p0|) 
and (|2T|) to search for EOSs yielding particular values of M^f . Having de- 



termined the k and 7 parameters for those EOSs yielding a particular M™™ , 
we then searched amongst those EOSs for the one which yielded the rotating 
star with the minimum period P m in- Figure ^| shows the results of our search 
for M^ x = 1.442 M . Each point along the curve defines an EOS which 
yields a maximum mass spherical star of mass M™% x = 1.442 M®. As one 
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moves along the curve towards smaller and smaller values of k, the minimum 
period allowed by the EOS continues to decrease and tends asymptotically to 
a minimum as k — > 0. 

That the period tends only asymptotically to a minimum as k — > can be 
understood by considering the EOSs sampled as one moves along the curve in 
Figure || The bottom plot of Figure || shows the EOSs sampled as one moves 
along the curve in Figure |6]. As one moves to small k the EOS is essentially 
unchanged, approaching the minimum period EOS (0) with constant pressure 
from p m to pc- In the limit as k one obtains the minimum period EOS 
(0) and the minimum period is given by the upper curve in Figure [| This 
is true for all reasonable M™^. 

5 DEPENDENCE ON MATCHING 
DENSITY 

Our results are fairly insensitive to the matching density at which the min- 
imum-period EOS d||) is matched to the FPS EOS for low densities. This is 
because the matching density is low compared to the central densities of the 
minimum-period stars, and as noted above, the causality limited region of 
the EOS determines the bulk properties of the minimum-period stars. Thus 
requiring that the star have a crust of normal matter and satisfy assump- 
tion (2) of section ^ increases the minimum period by less than 2.5% for a 
matching density of n m = 0.1 fm~ 3 and any . We further increased 

the matching density to n m = 0.25 fm~ 3 (well beyond where one trusts the 
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FPS EOS) and the minimum period again increased by less than 4.5 % for 
any M^f . This is shown in Figure [7]. Our current understanding of the low 
density EOS has a small effect on the minimum period determined only by 
the causality constraint (1) of section [2]. 

6 ACCURACY CHECK 

The numerical code used to construct rapidly rotating, fully relativistic com- 
pact star models was checked for accuracy in an extensive comparison to 
other codes (Eriguchi et al. 1996). The agreement of the codes for physi- 
cal quantities of specific models was 0.1 % to 0.01 % or better, depending 
on the stiffness of the EOS. Determining the maximum-mass model for a 
given EOS requires computing a large number of models until one computes 
a model within some tolerance of the true maximum mass and Kepler limit. 
Since the mass vs. central energy density curve is approximately flat near 
the maximum mass, we used a finely spaced grid and small tolerances (hence 
a large number of computed models) to determine the maximum-mass model 
for a given EOS. 

Still, for a given EOS the maximum-mass model is in general distinct 
from the minimum-period model. The intersection of the Kepler limiting 
curve with the line of onset of axisymmetric instabilities in a mass vs. cen- 
tral energy density plot defines the minimum-period model (Stergioulas & 
Friedman 1995; Cook et al. 1994). The onset of axisymmetric instabilities is 
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denned by the relation 

b =0 ' (22) 

where J is angular momentum. The central energy density at the intersec- 
tion may be higher or lower than the central energy density of the maximum- 
mass model. In general, however, the central energy density, mass, period 
and other characteristics of the minimum-period and maximum-mass models 
nearly coincide. We have confirmed that this is also true for the minimum- 
period EOSs ([l|) and @. Figure [8] is a plot of mass vs. central energy den- 
sity for sequences of models with constant angular momentum constructed 
using the minimum-period EOS (|]) and yielding M™% x = 1.442M . Also 
shown is the Kepler limiting curve and the axisymmetric instability line. 
The maximum-mass model has a central energy density of 4.9 x 10 15 g/cm 3 
and a period of 0.288 ms. The star having the maximum angular velocity for 
this equation of state has a central energy density of 5.3 x 10 15 g/cm 3 and a pe- 
riod of 0.282 ms, which is only 2% less than the period of the maximum-mass 
model. Computing the period of the maximum-mass model, rather than the 
true minimum- period model, saves a tremendous amount of computing time 
and introduces an error not larger than 2 %. 

7 COMPARISON TO EARLIER WORK 

Glendenning (1992) first estimated a causally limited period for gravitation- 
ally bound stars. Rather than numerically computing rapidly rotating mod- 
els, he computed nonrotating models and estimated the minimum period for 
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a given EOS using the empirical formula 



P ( M max \ ~^ / R ma:E \ f 

_= = 0.873 (-^) (^) . (23, 

This empirical formula involves only the mass M^f and radius of the 
maximum mass nonrotating star for the given EOS. Haensel and Zdunik 
(1989) and Friedman et al (Friedman, Ipser, & Parker 1989) constructed 
this formula using realistic EOSs; it has has an uncertainty of about 10 % 
(the Haensel-Zdunik coefficient is more accurate and has in these units the 
value 0.83). For a 1.442 M & maximum-mass nonrotating star Glendenning 
estimated the absolute minimum period to be 0.32-0.33 ms, which is 10-14 
% larger than our computed value of 0.28 ms. The difference is as large as 
the uncertainty of the empirical formula fl23|) , which is not surprising since 



(23) was constructed for a a set of EOSs that is vastly different from the 
minimum-period EOSs. Our numerical results yield an empirical formula for 
the class of minimum-period EOSs which is fairly insensitive on the matching 
number density n m . For a matching density of n m = 0.1/m~ 3 the empirical 
formula is 



P ( M max \ ~5 / Umax \ § 

rmin = o 740 sph ph (24) 

ms \ M & J \ 10 km) ' v ; 

with an approximate uncertainty of 0.3 %. This formula is accurate to about 
1% for 0.0 < n m < 0.25 fm~ 3 . 

Rather than identifying a single class of minimum-period EOSs, Glen- 
denning gave specific examples of minimum-period EOSs with and with- 
out a constant pressure region. As we have demonstrated, there is only 
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one class of minimum-period EOSs, and each EOS in the class has a con- 
stant pressure region. We believe that Glendenning missed this and pro- 
posed two classes of minimum period EOSs because he used a modified 
Levenberg-Marquardt method (Press et al. 1992) to minimize the function 
f(M, P) = wi(M — M^ff) 2 + w 2 (Pmin) 2 , where W\ and w 2 are weights. Our 
own trials with this method showed that it is difficult simultaneously to ob- 
tain high accuracy for both M™™ and P m i n - 

8 ABSOLUTE MAXIMUM MASS OF 
UNIFORMLY ROTATING NEUTRON 
STARS 

Among EOS satisfying assumptions (0-2) of section |2], that which yields a 
spherical model of maximum mass has a simple and unique form, consist- 
ing of two parts. The first part for low densities is a known low density 
EOS (earlier authors used the BPS (Baym, Pethick, &; Sutherland 1971) and 
Negele-Vautherin EOS). The second part, for densities greater than some 
matching density e m , is the causal limit EOS. Rhoades and Ruffini (1974) 
chose a particular e m and computed a maximum mass of 3.2 M Q . Although 
3.2 Mq continues to be quoted as the theoretical maximum mass for neutron 
stars, M max is sensitive to the matching energy density e m (as Hartle and 
Sabbadini (1977) pointed out) and a more accurate statement of the upper 
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limit set by causality on the mass of compact spherical stars is (Hartle 1978) 



This is in sharp contrast to the class of minimum-period EOSs, where due 
to the constant pressure region, the minimum period is insensitive to the 
matching density. So 3.2 M corresponds to a specific choice of e m and 
is not the theoretical maximum mass of neutron stars, which one needs to 
recompute whenever new certainty of the EOS allows one to increase the 
matching density. 

Friedman and Ipser (1992) assumed that the EOS limited by causality 
yields the maximum mass for rotating stars as it does for nonrotating stars. 
They computed, using an independent code, the maximum mass for different 
matching densities, assuming the BPS or Negele-Vautherin EOS for low den- 
sities. We also computed the maximum mass of rotating compact stars using 
the FPS EOS and our numerical code, which is based on a different algorithm 
than the Friedman and Ipser code. Our results are summarized as 



which is 3 % larger than the Friedman and Ipser result, and is within the 
numerical accuracy of the code used by Friedman and Ipser. Using the FPS 
rather than the BPS EOS did not significantly affect the maximum mass. 



the FPS EOS. If future research establishes the accuracy of FPS for higher 




(25) 




(26) 



In (|26|), 2 x 10 14 g/cm 3 is roughly the energy density up to which we trust 
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energy densities, then ( PS| ) can be used to give an updated, lower value for 

]\/fmax 
lvl rot ■ 
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Figure 1: Schematic representations of the minimum-period EOSs (|T]) and 
(0). The minimum-period EOS ([[]) does not match a known low density EOS; 
the pressure vanishes for e < ec, and is at the causal limit with dp/de = 1 for 
e > ec- The minimum-period EOS @ matches the FPS EOS to a constant 
pressure region at an energy density e TO . For e > ec the EOS is at the causal 
limit with dp/de = 1. Both axes are linear. 
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Figure 2: Minimum period EOSs ([I]) and (^) yielding a maximum-mass non- 
rotating star of M™^ x — 1-442 M & . The solid curve is the minimum period 
EOS ([[]), which does not match at low density to a known EOS. The energy 
density at which the EOS becomes causal is ec = 2.156 x 10 15 g/cm 3 . The cen- 
tral energy density of the minimum-period star is 4.778 x 10 15 g/cm 3 and is in- 
dicated by a filled circle on the solid curve. The dotted curve is the minimum- 
period EOS (H]), which matches to the FPS EOS at low density. The energy 
density at which the EOS becomes causal is 6q = 2.157 x 10 l5 g/cm 3 . Above 
this density the EOSs are nearly identical and the solid curve overlaps the 
dashed curve. The central energy density of the minimum-period star is 
5.274 x 10 15 g/cm 3 and is indicated by a filled triangle on the dotted curve. 
The matching number density is n m = 0.1 fm~ 3 . Both axes are logarithmic. 

Figure 3: The minimum period P m in allowed for a rotating, relativistic star 
as a function of the mass M™^ of the maximum-mass spherical star allowed 
by the EOS of the stellar matter. The upper curve is constructed using the 
minimum-period EOS @, which matches at low density to the FPS EOS. 
The lower curve is constructed using the minimum-period EOS ([J), which 
does not match at low density to a known low density EOS. Each curve 
divides the mass-period plane into two regions. The region below a curve 
is not accessible to stars with EOSs that satisfy the assumptions listed in 
section |2], and hence the region below the lower curve is not accessible to 
stars composed of matter obeying causality. 
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Figure 4: Minimum period P m in of rotating stars as a function of the allowed 
maximum mass M™£ x of spherical stars for EOSs which are perturbations to 
the minimum-period EOS @. The squares, triangles, and crosses represent 
EOSs obtained by different perturbations to the minimum period EOS, which 
are detailed in section |4.1| . The solid curve represents the minimum-period 
EOS (fj) and is the same as the upper curve in Figure fj. All of the points 
obtained from perturbed minimum-period EOSs lie above the solid curve. 



Figure 5: Both plots show EOSs for different k and 7 using the param- 
eterization (|2~0D and (|2l|) . The K'j parameterized EOS is fairly robust and 
can produce both stiff and soft EOSs. The bottom plot shows a sequence of 
EOSs, all of which yield M™H* = 1.442M . The period of the fastest rotating 
star decreases from left to right as the EOS approaches the minimum-period 
EOS (0). All axes are logarithmic. 



Figure 6: Minimum period allowed by EOSs given by ( |20|) and (|21|) for dif- 
ferent values of the parameter k. For given k, the parameter 7 is determined 
so that all EOSs sampled as one moves along the curve yield a maximum- 
mass nonrotating star with M™ff = 1.442 M . As k — > the period tends 
asymptotically to its minimum value. 
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Figure 7: The minimum period P m in as a function of the mass M™jf of the 
maximum- mass spherical star, for different matching number densities n m . 
From top to bottom the curves represent n m = 0.5,0.25,0.1, and 0.0 fm~ 3 . 
Note that for n m = 0.5 fm~ 3 the maximum- mass spherical model cannot 
exceed 2.38 M G) . 



Figure 8: Mass vs. central energy density for sequences of models with 
constant angular momentum, constructed using the minimum-period EOS 
([]]) with M™ff = 1.442 M Q . The solid curves are sequences of models with 
constant angular momentum, the bottom one corresponding to nonrotating 
models. The dashed curve is the sequence of models rotating at the Kepler 
limit while the dotted curve is the axisymmetric instability line defined by 



(|22f) . The maximum- mass model is distinct from the model with maximum 
angular velocity VL max . The period of the maximum-mass model is only 2% 
less than the period of the fastest rotating model. 
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